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Abstract 

We extract values for the free symmetry energy as a function of the fragment size (the proton 
number Z) from antisymmetrized molecular dynamics (AMD) calculations of calcium collisions. 
Simple statistical physics describe well the distribution of hot nuclei at breakup, provided the 
surface symmetry term in the free energy is much smaller at high excitation than in ground state 
nuclei. This result may reflect the condition of low density and finite temperature when these 
systems disassemble. 

PACS numbers: 25.70.Pq 
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Free energies play an important role in mixed phase environments. For nucleonic systems, 
such environments may be found in subsaturation density systems formed in nucleus-nucleus 
collisions, supernova collapses ll| and in the inner crusts of neutron stars In these 

scenarios, the free energy of fragments largely defines the balance between denser fragments 
and the more dilute nucleonic gas. Key uncertainties to the prediction of free energies are 
their dependencies on temperature, density and isospin asymmetry. The latter uncertainty 
is particularly relevant to the inner crust of neutron stars, where the nuclear isospin asym- 
metry term stabilizes nuclear droplets enveloped by a lower density neutron gas 0,0]. The 
development of experiment constraints on the asymmetry term constitutes an important 
scientific objective. 

Here we focus on multifragmentation in nucleus-nucleus collisions. Many aspects of such 
data have bee. .uecessMly desenbed bv e,u„ib.„. .a«s*a> .node. *at a.„.e 

fragments to be produced by a low density phase transition that occurs during the expansion 
stage of a collision. The success of such descriptions suggests applications of statistical 
theory to such collisions may provide constraints on the free energies and their asymmetry 
dependence. Such comparisons provide an experimental test of the validity of equilibrium 
descriptions. One can also test the assumption of local thermal equilibrium by comparing 
the predictions of modern transport theories to equilibrium models. Calculations of isotopic 
observables are particularly sensitive to the degree of chemical equilibrium achieved. 

We investigate these issues with the antisymmetrized molecular dynamics (AMD) code 
of Ref. 0. In the present work, we analyse the yields of fragments produced in AMD 
simulations for central collisions of nuclei with 40 < A < 60 at an incident energy E/A = 35 
MeV. Even though AMD does not require the assumption of equilibrium, we have shown in 
a previous paper P] that some observables obtained from predicted isotope yields, such as 
isoscaling, are consistent with the expectations from a statistical interpretation of fragment 
formation. 

If the yield of fragments, produced with N neutrons and Z protons in these collisions, is 
governed by a statistical process at constant pressure, the fragment yield Y{N, Z) can be 
related to the nuclear free energy G^udN, Z) by the relation 



F(iV, Z) ocexp 



(1) 



T T T 

where fin and /ip are the neutron and proton chemical potentials. At temperatures relevant 
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FIG. 1: Binding energies (not per nucleon) of nuclei calculated by AMD with the Gogny force. 
The quantities — BE(A^, Z) + 8(A^ + Z) MeV are shown by filled and open squares for even-Z and 
odd-Z nuclei, respectively. Lines connect isotopes. 

for multi-fragmentation, the temperature dependence of the asymmetry term in the free 
energy for bulk nuclear matter is expected to be negligible. However, the temperature 
dependence of the surface symmetry term of the free energy is not well known. The present 
study suggests that the surface contribution to the symmetry free energy may be weakened 
under the conditions associated with fragment formation. 

Calculations within the AMD formalism have demonstrated its capability of reproducing 
many quantal features of ground state nuclei . The present AMD code was streamlined for 
collisions; it nevetheless reproduces ground state binding energies of nuclei to within about 
0.5 MeV/nucleon. Figure ^ shows the binding energies of nuclei BE(A^, Z) with A < AO 
which are obtained by minimizing the energy within AMD by adopting the Gogny force 
jlfll |. In order to remove the strong A-dependence, we plot — BE(A^, Z) + 8A MeV. We 
note that the free energy equals the energy itself at zero temperature and pressure, i.e. 
Gnuc{N, Z,T = 0,P = 0) = -BE{N, Z). 

Much of our present knowledge of the symmetry energy has been extracted by fitting 
the nuclear binding energies with the liquid-drop mass formula which contains a symmetry 
energy term such as 

E^,^{N,Z) = c{A) ^^~^^^\ (2) 

In the simplest mass formula the coefficient is independent of the nuclear size, c{A) = 
Csym, which assumes the volume nature of the symmetry energy. On the other hand, advanced 



mass formulas 



have introduced the A-dependence of c{A) as the surface effect. 
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The extraction of c{A) from the nuclear binding energies is not quite straightforward even for 
ground state nuclei. For example, if one extracts the symmetry energy by using the energy 
difference of neighboring nuclei |l5|, the symmetry energy is largely a fluctuating function 
in the nuclear chart due to shell and paring effects. Nevertheless, the global fitting with the 
assumption of c{A) = + CsA'^^^, together with the standard volume, surface. Coulomb 
and paring terms, usually results in a reasonable value of the coefficients. For example, we 
obtain Cv = 27.3 MeV and Cs = —23.7 MeV by fitting all the available information of the 



binding energies of nuclei [16|, while we obtain = 30.1 MeV and Cg = —35.1 MeV by 
fitting the nuclei with 7 < A < AO. 

If we fit the AMD binding energies in Fig. ^for 7 < A < 40, we obtain = 30.9 MeV 
and Cs = —35.2 MeV. The extracted value of Cy is comparable with the symmetry energy 
in infinite nuclear matter at saturation density po- To a good approximation, the nuclear 
matter EOS at zero temperature can be written as 

E{p, 5) /A = Eip, S = 0)/A + C3ym(p)5', (3) 

where p = Pn + Pp and 6 = {pn — Pp) / p, and the symmetry energy coefficient has the value 
C'sym(po) = 30.7 MeV for the Gogny force. The present study of heavy ion fragmentation 
reactions will provide information about the symmetry free energy at subsaturation densities 
and finite temperatures. 

In order to extract the symmetry free energy from Eq. (^, it is necessary to obtain 
the fragment yields Y{N, Z) for a wide range of isotopes. This is practically impossible if 
we utilize the products produced in only one reaction system. However, we overcome this 
limitation by combining the results from various reaction systems here labeled by an index 
I. We calculate ^oCa + ^OQa {i = l), ^^Ca + ^^Ca (i = 2), ^^Ca + ^^Ca {i = 3) and ^^Pe + ^^pg 
{i = 4) collisions at zero impact parameter and an incident energy E/A = 35 MeV. Previous 
calculations with an equivalent version of AMD and the Gogny force [itI . Iisl reproduce 
experimental data for various fragment observables in ^°Ca + ^°Ca collisions at E/A = 35 
MeV. 

AMD represents the wavefunction of the colliding system by fully antisymmetrized prod- 
ucts of Gaussian nucleon wave packets and propogates these wave packets microscopically 
during the collision [l, 12, 3|- The centroids of the nucleonic wave packets move determin- 



istically through the mean field potential formed by the interactions with other nucleons. In 
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addition, the followed state of the simulation branches stochastically and successively into 
a huge number of reaction channels. The branching is caused by the two-nucleon collisions 
and by the splittings of the wave packets. The interactions are parameterized in the AMD 
model in terms of an effective force acting between nucleons and in terms of the two-nucleon 
co.„.o„ CO. secfo... We adopt the Oo,.y effective fo.ee Q in th. pape. 

We simulate collisions by boosting two nuclei whose centers were separated by 9 fm and 
calculating the dynamical evolution of each collision until t = 300 fm/c. The numbers of 
simulated events are 1040, 949, 978 and 1400, respectively, for the four systems. In central 
collisions, as shown in previous papers the projectile and target basically penetrate 

each other and many fragments are formed not only from the projectilelike and targetlike 
parts but also from a neck region between the two parts. The nuclear matter seems to be 
strongly expanding one-dimensionally in the beam direction. 

As we have shown in Ref. j^, the AMD results satisfy the isoscaling relation 

F,(Ar,Z)/r,,(iV,Z)oce-^+^^ (4) 

with parameters a and f3, for the fragment yields from two different reaction systems {i 
and i') with different neutron-to-protion ratios. When we have four reaction systems {i = 
1,2,3,4), the fragment yields Yi{N,Z) are related to the yields Yi{N,Z) of the reference 
system by 

ri(iV, Z) = Yi{N, Z)e-''^^-^'^^\ (5) 

The yields Yi {N, Z) for one reaction system can have good statistics only in a narrow region 
of (A^, Z). The isoscaling relation is utilized with respect to A^ to extend the range of (A^, Z). 
We take the ^°Ca + ^°Ca system as the reference system and therefore ai = 'ji{Z) = 0. The 
quantity — In Yi(A^, Z), i.e. the fragment yields for ^°Ca + ^°Ca, can be represented in four 
different ways by using the scaled yields for different reaction systems which have good 
statistics in different regions. They are combined by 

4 

K{N, Z) = J2 Z)[- lnY,{N, Z) + a,N + 7,(Z)] (6) 

i=l 

where the averaging weights Wi{N,Z) are determined by minimizing the statistical errors 
in K{N, Z) for individual (A^, Z). The isoscaling parameters a,, which is common to all Z, 
have been obtained by the isoscaling fit The parameters 'jii^Z) {i = 2,3,4) for each Z 
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FIG. 2: The values of K{N, Z) for 3 < Z < 18 are shown by symbols for the abscissa of + Z. 
The values are obtained by combining the results of '^'^Ca + ^'^Ca, "^^Ca + '^'^Ca, ^'^Ca + ^'^Ca and 
^^Fe + ^^Fe simulations. The error bars show the statistical uncertainty due to the finite number 
of events. The curve for each Z was obtained by fitting K{N, Z) using Eq. 0. 

are determined by optimizing the agreement of the quantities [—\iaYi{N, Z) + + 7i(2')] 
from different reactions i = 1,2,3 and 4. 

The values of K{N, Z) obtained from the AMD simulations are shown by the open and 
filled squares in Fig. |21 for odd- and even-Z nuclei, respectively. Even though the event 
numbers are not very large, K{N, Z) has been obtained for a wide region of (A^, Z). In most 
cases K{N, Z) covers more than 10 isotopes with good statistical precisions for each Z. 

So far, we have not assumed any specific form of K{N, Z). Nevertheless, the results show 
a very smooth behavior of K{N, Z) as a function of and Z. The shell and paring effects 
are weak in K{N, Z) compared to the ground state binding energies shown in Fig. ^ Each 
curve in Fig. |21 shows the fitting of K{N, Z) for each Z by a function 



where ^(Z), 1]{Z) and ({Z) are the fitting parameters. The result of the simulations is fitted 
well by this functional form. We choose the quadratic term similar to Eq. Q for convenience, 
so that the parameter C{Z) is directly related to the symmetry energy as shown below. 

The obtained values of C{Z) are shown in Fig.Olby solid points. Except for light fragments 
{Z < 5), C{Z) is a smooth function of Z which depends on Z very weakly. Within the 
uncertainty in the fitting procedure, the trend does not depend on the number of reaction 
systems included in Fig. |21 As expected, the Z-dependence becomes more smooth and stable 



K{N,Z) = aZ)N + r]{Z) + aZ) 



{N - Zf 
N + Z 



(7) 
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FIG. 3: The solid points are the extracted values of the coefficients C(^) of Eq. ((7| using the 
combined fragment yields of four systems shown in Fig. |2 The thick solid curve, the thin sold 
curve and the dashed line show functions C{^) oc 1 — k{2Z)^^/^ normalized at Z = 10 for k = 1.14, 
0.5 and 0, respectively. 
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FIG. 4: The solid points show C{Z)T when T = 3.4 MeV is assumed. The dot-dashed horizontal 
line shows the volume symmetry energy Cv = 30.9 MeV for the ground state nuclei calculated with 
AMD. The solid line shows the symmetry energy c{A = 2Z) = Cv + Cs(2Z)-^/3 for the ground state 
nuclei (cv = 30.9 MeV and Cs = -35.2 MeV). 
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by including more reaction systems (up to four in the present work) due to the extension of 
the range of A^. 

Next we will show that the extracted value of C{Z) can be regarded as the ratio of the 
symmetry energy to the temperature. Under the assumption of equilibrium, K{N, Z) is 
identified with the exponent of Eq. ^ with inverted signs. If we approximate the nuclear 
free energy G'nuc(^7 Z) by a smooth function of (A^, Z), K{N, Z) can be written as 

K{N, Z) = f{A, Z) + g{A)^^^ + aN + bZ, (8) 

where f{A, Z) and g{A) have smooth A dependence. Specifically, we associate the symmetry 
energy c{A) in Grmc{N^ Z) with g{A) by g{A) = c{A)/T. f{A, Z) and the linear terms in 
and Z can be associated with terms such as volume, surface and Coulomb terms in the free 
energy and the chemical potential terms in Eq. (^. Because the range of important for 
each Z is limited, we can expand f{A, Z) and g{A) with respect to A^ for each Z. By denoting 
the typical value of A for each Z by A{Z), f{A, Z) ^ f{A{Z), Z)+f'{A{Z), Z){N+Z-A{Z)) 
and g{A) ^ g{A{Z)). The higher order terms are negligible as can be checked directly by 
assuming the nominal values of the liquid-drop coefficients. 

Thus we have shown that Eq. (|7j) is a natural form and that the coefficient C{Z) is related 
to the symmetry free energy c{A) by 

aZ) = c{A{Z))/T. (9) 

In principle, ({Z) could depend on Z, reflecting the size-dependence of the symmetry energy 
c{A). However, the extracted values of C{Z) in Fig.Elare almost independent of Z for Z > 5. 
This suggests a reduction of the surface contribution to the symmetry energy. The three 
curves in Fig. IHlshow functions CiZ) oc 1 — k(2Z)^^^^ for different surface-to-volume ratios 
A; = 1.14 (thick solid line), 0.5 (thin solid line) and (horizontal dashed line), respectively. 
The curves are normalized a.t Z = 10. The Z-dependence of C{Z) cannot be explained by 
the surface-to- volume ratio k = — Cg/cv = 1.14 for the symmetry energy of ground state 
nuclei. The result shows that the surface effect in ({Z) is reduced to between k = and 0.5. 

There can be several possible explanations for the weakening of the surface symmetry 
free energy. First of all, it is not very surprising that the coefficient in the free energy at 
a finite temperature is different from that at the zero temperature. As is well known, the 
surface tension reduces towards zero when the temperature is raised towards the critical 



8 



temperature. A similar effect 



Fermi surface calculations 
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las been obtained for the symmetry free energy in Thomas 



221 . If we adopt the formula in Ref . 



2jJ , the reduction factor 



of the surface symmetry energy is approximately [1 — (T/Tc)^]^ = 0.91 for T = 3.4 MeV (see 
below) and the critical temperature = 16 MeV. This reduction is, however, not sufficient 
to explain the weak Z-dependence of C{Z). Another explanation may be associated with the 
fact that fragments are not isolated when they are formed. When the density fluctuation 
is developing from a uniform low density matter, the fragments are still interacting with 
attractive force through their surfaces. Therefore, surface free energies could be expected 
to be smaller for these fragments than for totally isolated fragments. Independent of the 
physical origin for the weakening of the surface symmetry free energy, it suggests that the 
volume quantity, which is the same as that in the inflnite nuclear matter, can be directly 
obtained by the analysis of the fragmentation results even though the produced fragments 
are not very large. 

In the above analysis, only the ratio of the symmetry free energy to the temperature 
is obtained. In the present theoretical approach, it is possible to get the density and the 
temperature by studying the response of the results to a change of the symmetry term in the 
effective force. In Ref. we have derived the density p ~ 0.08 fm"'^ and the temperature 
T ~ 3.4 MeV for the same reaction systems studied in the present paper. This result is 
consistent to the idea that the fragments are formed in a low density nuclear matter at a 
finite temperature. The solid points in Fig. EJshow C{Z)T so that the extracted C{Z) can 
be directly compared to the symmetry energy. The value of the symmetry energy extracted 
from fragmentation is significantly lower than that of the bulk nuclear matter at normal 
nuclear matter density, Csym(po) ~ Cv = 30.9 MeV shown by the dot-dashed horizontal fine 
in Fig. m and this must be due to the reduced density, p ~ 0.08 fm"^ < pq. 

The solid line in Fig. |3] shows the symmetry energy [30.9 - 35.2(2Z)-i/3] ^eV that 
has been obtained by fitting the ground state binding energies calculated with AMD. The 
extracted symmetry energy is smaller than the ground state symmetry energy probably 
because of the greater reduction from the surface contribution when fragments are formed. 

In conclusion, we have analysed AMD simulation results from nuclear collisions of vari- 
ous nuclei with different neutron-to-proton ratios. Fragment yields from different reaction 
systems are combined using the isoscaling relation. The availability of fragment yields over 
a wide range of N and Z allows us to extract the symmetry energy at low density when 
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fragments are formed. The results are consistent to the idea that the fragments are formed 
in nuclear matter at a low density and at a finite temperature. The extracted symmetry 
energy shows almost no surface effect in it, which suggests that the property of infinite 
nuclear matter can be directly obtained from the information of fragmentation. 
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